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Strained Fe-Co films have recently been demonstrated 
to exhibit a large magnetocrystalline anisotropy (MCA) 
and thus to be of potential interest as magnetic storage 
material. Here, we show by means of density-functional 
(DF) calculations, that chemical order can remarkably 
enhance the MCA. 



We also investigate the effect of relaxation perpendicu- 
lar to the applied strain and evaluate the strain energy 
as a function of Co concentration and substrate lattice 
parameter. On this basis, favourable preparation routes 
for films with a large perpendicular anisotropy are sug- 
gested. 
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1 Introduction Materials for high-density magnetic 
recording media (hard disc drives) have to obey two com- 
peting requirements (T\. On the one hand, stability against 
thermally activated switching is guaranteed only if the 
volume-integrated MCA of a single storage bit is larger 
than about 50 fceT, where fee is the Boltzmann constant 
and T is the operating temperature Ul- On the other hand, 
the magnetic field needed to write a bit is, apart from de- 
magnetisation effects, proportional to K^/Ms 13], where 
Ku denotes the (uniaxial) MCA energy and Ms is the sat- 
uration magnetisation. Since the write field is constrained 
by the construction of the write head a large of 
the storage material is desirable along with a large Ku, the 
latter demand resulting from the stability requirement. 

Fe-Co alloys are well-known to have a high saturation 
magnetisation at ambient conditions, but their bulk MCA 
is minute due to their cubic symmetry. It is however pos- 
sible to achieve high anisotropy values without sacrificing 
the advantageously large magnetisation by designing arti- 
ficial structures that break the cubic symmetry. One way 
to achieve this goal would be to manufacture Fe|Co su- 
perlattices. It has been demonstrated experimentally and 
confirmed by DF calculations, that [(1 10)-Fe|Co]„ super- 
lattices show a large in-plane magnetic anisotropy 0. Un- 
fortunately, advanced recording techniques require a mag- 



netisation orientation perpendicular to the surface, i.e., an 
out-of-plane anisotropy. A second possibility would be the 
preparation of ultrathin films. Monolayers of Fe-Co on 
Pt(lll) show a maximum out-of-plane anisotropy energy 
of about 500 /zeV per atom both in DF theory and in low 
temperature experiments. However, they have Curie tem- 
peratures close to room temperature [0 and are, thus, not 
suited for the discussed application. 

A third, recently suggested route is to fabricate 
strained, bulk-like Fei_^Coj; films by epitaxial growth 
on a suitable substrate [:6:|. This idea relies on the fact 
that alloys which are cubic in their bulk phase can be 
grown as metastable tetragonal films, if they are deposited 
on substrates with a fourfold symmetric surface, e.g., the 
(001) surface of cubic crystals. It was shown in a recent 
study that the Feo.yPdo.a alloy can be prepared to form 
50 nm thick, i.e., bulk-like, epitaxial films on a number of 
substrates. Thus, the centred-tetragonal structure of these 
films with lattice-parameter ratios c/a between 1.09 and 
1.39 spans almost the whole range from BCC (c/a = 1) to 
FCC (c/a = s/2) Q. 

Turning back to Fei-^^Co^; alloys, the out-of-plane val- 
ues of AEmca = Ku {T = 0) predicted by DF calcu- 
lations reach up to 800 /ieV per atom for x ~ 0.6 and 
c/a between 1.20 and 1.25 IS). We will demonstrate below, 
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that the elastic strain energy of Fei-^^Coa. is small enough 
in the interesting region of x to allow the stabilisation of 
films with a decent thickness. (Note, that quantum-size ef- 
fects strongly influence the MCA of thin films at least up 
to 10 monolayers (ML) or about 1.5 nm thickness fS).) 

Meanwhile, the mentioned prediction by Burkert et al. 
(abbreviated below BNEH 16]) has been confirmed by sev- 
eral experiments ||9][T0][TT][T2] ■ These confirmations should 
be viewed with caution, since in all quantified cases the 
measured values of Ku seem to be lower than the pre- 
dictions for comparable parameter values of c/a and of 
X by factors of two to four In detail, A'„ = 108 /xeV 
per atom was found for Feo.5Coo.5|Pd(001) films with 3- 
10 ML thickness, about half the predicted value llOl ; for 
[Feo.36Coo.64|Pt(001)],i superlattices, only one half ||9l or 
one quarter rtlll of the theoretical value was measured, de- 
pending on the way of comparison (dedicated calculation 
for a certain superlattice f9] or evaluation of bulk data from 
various superlattice data 111 11 ). One should note that in the 
latter case the anisotropy was measured at room tempera- 
ture, which explains a part of the discrepancy. 

Looking for a reason of the possibly systematic dis- 
agreement between theory and experiment, a first thought 
might blame the mentioned quantum oscillations. How- 
ever, data evaluation in two of the three mentioned cases 
considered fits to a number of systems with varying film 
thickness, and the third comparison involved a calculation 
for the specific geometry. Another possible source of de- 
viation might be the local density approximation (LSDA) 
applied by BNEH. We argue, that this is unlikely since in 
most known cases LSDA results for AEmca underesti- 
mate the experimental i^Tu-values. 

Here, we are going to advocate a third idea: While 
the experiments in all cases known to us were performed 
on presumably chemically disordered (though structurally 
well-ordered) alloys, the calculations by BNEH employed 
the so-called virtual crystal approximation (VCA) to simu- 
late the alloy. In this approximation, the electron number is 
adjusted to its correct value according to the alloy compo- 
sition by allowing the atomic nuclei to carry a non-integer 
charge. In this way, a perfect chemical order is introduced. 
In our particular case, all Fe and Co atoms are replaced by 
only one kind of atoms with atomic charge 26 + x. This 
means, the VCA describes a chemically ordered structure 
with the correct electron number of the chemically disor- 
dered alloy. 

As the main result of the present work, we will show 
that strained Fe-Co films with chemical order can have 
a much larger magnetic anisotropy than chemically dis- 
ordered films. Beyond this main point, we will contrast 
the strain-energy landscape with the related magnetocrys- 
talline anisotropy energy in order to find promising prepa- 
ration parameters for films with a large perpendicular 
anisotropy. 

The following sections contain computational details, 
numerical results and their discussion, and the conclusions. 
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Figure 1 Geometries used in the calculations: Constant 
volume Bain path using the experimental volume of BCC 
iron (thick dashed lines. Model I) and epitaxial Bain path 
calculated in VCA for different cobalt concentrations x 
(symbols. Model II). The shaded area with dashed mar- 
gins for .T = and x ~ 1 denotes an a;-weighted aver- 
age between the calculated volumes of BCC iron and FCC 
cobalt. Thick full lines indicate c/a ratios of -\/2 (FCC) and 
1 (BCC). 

2 Computational details The DF calculations were 
carried out with the all-electron, full-potential local-orbital 
code FPLO, version 8.00-31 11311141 . The generalised gra- 
dient approximation (GGA) in the parameterisation by 
Perdew, Burke, and Ernzerhof iflSl was used. We preferred 
this approach against the LSDA previously used by BNEH, 
since we are interested in the evaluation of strain energies 
in addition to the magnetic quantities. The valence basis set 
comprised 3s, 3p, id, 4s, Ap, Ad, and 5s states. Brillouin 
zone integrations were performed using the linear tetrahe- 
dron method with Blochl corrections. Structural properties 
and spin magnetism were evaluated in a scalar relativistic 
mode, while a four-component fully relativistic mode was 
used for the calculation of the MCA. 

Five different structural models were employed to de- 
scribe the atomistic geometry and the chemical order of 
the strained alloy films. All these models are based on the 
assumptions that (i) the films are thick enough to disregard 
the influence of the surface and of the substrate on the elec- 
tronic structure and that (ii) the substrate-film interaction is 
nonetheless strong enough to let the in-plane lattice param- 
eter be dictated by the substrate. These assumptions are in- 
herent to the concept of the epitaxial Bain path (EBP) lfT6l 
[TTl that we applied in a part of the calculations. 

Model I: A centred-tetragonal (BCT) structure with 
one kind of atoms of atomic number 26 + x (VCA). The 
lattice parameter c was determined by the assumption of 
constant volume, where the latter was chosen to be the ex- 
perimental volume of BCC Fe at room temperature, i.e., 
c = 2 X 11.78 A? /a? (Figure 1, thick dashed Une, the so- 
called constant volume Bain path, CVBP). 
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Model II: Same structure (BCT) and chemistry (VCA) 
as in Model I, but the lattice parameter c was calculated 
by minimisation of the total energy E{a, c, x) for a and x 
fixed, using spin polarised GGA. Related data are given in 
Figure 1 with symbols for different x. These data constitute 
the EBP. For comparison, CVBPs are given (dashed lines 
including shaded area) using the GGA volumes of BCC Fe 
{x = 0) and of FCC Co (x = 1). The EBP values of c/a are 
close to the related CVBP values except for small x in the 
region close to the FCC structure. This can be understood 
by the strong x-dependence of the spin moment /.is in this 
region, see Figure |4]below. 

Contour plots presented below are based on data cal- 
culated with Model I or II on grids with Ax = 0.1 and 
Aa = 0.05 A. 

Model III: LIq structure with stacks of single quadratic 
Fe- and Co-layers, x ~ 0.5. This structure is the strained 
variant of the known a'-Fe-Co (B2) phase, where c/a = 1. 
The lattice parameter c was calculated by minimisation of 
the total energy E{a, c) for a fixed (EBP). 

Model IV: Stacks of single FCC-like (Fe, Co)- and 
(Co, Co)-layers, x ~ 0.75, and EBP condition for c. This 
structure is derived from the LI 2 structure which is ob- 
tained for c/a = V2- 

Model V: A 2 X 2 X 2 BCT supercell with 16 atoms, 
used to describe chemical disorder by means of an ensem- 
ble average. FegCog {x = 0.5), FegCoio (x = 0.625), and 
Fe4Coi2 {x = 0.75) were considered. The atoms were ar- 
ranged such that the nearest neighbour patterns match the 
completely disordered alloy as closely as possible. Starting 
from a single atomic configuration, a symmetry-adapted 
ensemble average was constructed in such a way that its 
MCA correctly vanishes for a = c (BCC structure). De- 
tails of the method are described in Reference JTS], where 
32-atom supercells modelled the case of a slightly distorted 
FCC structure. 

The space groups 139, 139, 123, 123, 1 were used in 
the calculations for structure Models I ... V, respectively. 
Structure optimisation was carried out for Models I ... IV, 
where fc-meshes with 24 x 24 x 24 points in the full Bril- 
louin zone (BZ) for all models were used. For the more 
complex Model V, internal relaxation is expected to be of 
minor importance, since Fe and Co atoms have almost the 
same atomic volumes. Thus, we abstained from a relax- 
ation of the inner degrees of freedom. Further, the atomic 
volume obtained with Model II at the appropriate x-value 
and a = 2.65 A (a value of particular interest, see next 
section) was taken to construct a CVBP for Model V. 

For the Models I ... IV, the MCA energy was evaluated 
from independent self-consistent total energy calculations 
for magnetic moment orientations along the a;-axis (i?ioo) 
and along the z-axis (i^ooi). 

^-E^MCA = £^100 — £-001 ■ (1) 

The in-plane anisotropy (i?ioo — ^-iio) was not considered. 
In these calculations, finer /c-meshes than in the structure 



optimisations were used: 48 x 48 x 48 points in the full 
BZ for all models. 

To cope with Model V, we relied on the so-called 
magnetic force-theorem 1191 : the total energy difference 
^^'mca was approximated by the difference of band en- 
ergy sums, evaluated with 8x8x8 fc-points in the full 
BZ. In Eqn. (1) i?ioo is to be replaced by (i^ioo + ^'Oio)/2, 
where -Eqio is the band energy calculated for a magnetic 
moment orientation along the y-axis. 

3 Results and discussion 



E — Eq [meV/atom] 









5/ ; 1 ; 
















40 




iO + 






















i 

1 1 ', 


■ 1 







0.2 0.4 


0.6 


0.8 



X 



Figure 2 Calculated strain energy (Model II) as a func- 
tion of the cobalt concentration x and of the in-plane lat- 
tice parameter a. For each value of x, the total energy is 
referred to the lowest energy Eq{x) along the respective 
EBP, i.e., there is a line ao{x) for which the presented en- 
ergy is zero (not shown). The (unstrained) ground state is 
BCC for X < 0.76 and FCC for x > 0.76. (Due to interpo- 
lation from a grid the jump from BCC to FCC at 2; = 0.76 
slightly deviates from being vertical in the picture.) The 
symbol + labels the point of the largest MCA according to 
Fig. Eh. 

3.1 Strain energy Figure |2] shows a contour plot of 
the strain energy, 

EstTs,inia,x) ^ E{a,x)-mmE{a,x) := E{a, x)-Eo{x) , 

a 

(2) 

with x-dependent ground-state energy Eq, evaluated with 
Model II (VCA, EBP). There are two valleys correspond- 
ing to the BCC (x < 0.76) and FCC (x > 0.76) ground 
state, separated by a saddle point at x = 0.78; a = 2.66 A. 
The calculated critical Co concentration of about 76% for 
the FCC-BCC transition is consistent with the known ex- 
istence range of the BCC-like Fe-Co B2 phase up to 72% 
Co and two-phase behaviour (FCC + BCC) at higher Co 
concentration 1201 . 

The strain energy is an important quantity that deter- 
mines the feasibility of epitaxial growth. In a simplified 
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picture, the achievable thickness of a metastable film is in- 
versely proportional to iSstrain- Recent work on Feo.rPdo.s 
films demonstrated the possibility to grow 50 nm films 
(about 300 ML) with a strain energy of 6 meV per atom, 
estimated by DF calculations l?]. Considering Fe-Co films 
on a Rh(OOl) substrate (a = 2.69 A), we predict a strain 
energy between 40 and 80 meV per atom for x < 0.65. 
Indeed, experiments indicate that epitaxial growth is stable 
up to about 15 ML, while thicker films seem to become 
crystallographically disordered llTOl . 

From this perspective, it seems worthwhile to try the 
growth of thicker films in the saddle point region of Figure 
|2] The related in-plane lattice parameters define a region 
of very high strain (either BCC or FCC serve as zero-strain 
reference), but the strain energy is comparably low. As an 
additional advantage of this area the in-plane stress is low 
(zero at the saddle point) which might additionally facili- 
tate the epitaxial growth. We note, that in reality the saddle 
point area might appear at about 0.04 A higher a-values, 
since the GGA calculations find somewhat smaller ground- 
state volumes than known from experiment, see Figure [T] 



ure |2]l of 20 . . . 25 meV per atom. We conclude, that the 
a-dependence of the strain energy is mainly determined by 
the electron number featured in the VCA. Structural de- 
tails and chemical disorder only slightly modulate the be- 
haviour. 
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Figure 3 Calculated strain energy for cobalt concentra- 
tion X ~ 0.75 as a function of the in-plane lattice pa- 
rameter a. Results for two different chemically ordered 
structures (Model II, denoted VCA, and Model IV, denoted 
LI 2) are compared with results for ensemble-averaged su- 
percells simulating chemical disorder (Model V, denoted 
EAVG). The total energies are referred to the respective 
lowest value Eo{x — 0.75). The VCA data are the same as 
those used in Figure |2] but were evaluated here on a finer 
grid. 

To check, whether different types of structural or chem- 
ical order have an influence on the film stability, we com- 
pare strain energies obtained with three different struc- 
ture models, Figure|3] The calculations were performed for 
X = 0.75, a line cutting the a — x-plane close to the sad- 
dle point. All three models show very similar dependence 
of i?strain on fl, with minima at the FCC- and BCC-(like)- 
structures and a maximum (close to the saddle point in Fig- 



Figure 4 Calculated spin magnetic moment (in VCA, 
Model II) as a function of cobalt concentration x and in- 
plane lattice parameter a. The symbol + labels the point of 
the largest MCA according to Fig. |5j). 

3.2 Spin moment For magnetic storage applications, 
a large saturation magnetisation is needed to allow switch- 
ing with a magnetic field of limited strength. Figure |4] 
shows the VCA spin moment, which accounts for about 
90% im of Ms (T = 0), as a function of x and a. We 
find high values of /i,, depending only weakly on x and 
being almost insensitive to a variation of a, except in the 
region x < 0.2; a < 2.60 A. The weak x-dependence 
can be understood by assuming constant atomic moments 
for X > 0.2 and transition from weak to strong ferro- 
magnetism of iron in the region < a; < 0.2 (Slater- 
Pauling behaviour). For the case of Fe-Co films on Cu(OOl) 
(a = 2.55 A), this behaviour has been confirmed in exper- 
iment II22I . 

The pronounced moment reduction in the Fe-rich area 
below a = 2.60 A is related to the instability of fer- 
romagnetism in FCC iron. In the present GGA calcula- 
tions, restricted to collinear ferromagnetic states, we find 
a low-moment solution for very small Co concentrations 
and a high-moment solution for x w 0.2. This strong x- 
dependence of the spin moment has an effect on the EBP 
(Figure [Til by magneto-volume coupling: for a; = 0, the 
volume is reduced in comparison to the constant-volume 
assumption, while it is considerably enhanced for x ~ 0.2. 

The earlier results for fig (x, c/a) by BNEH agree qual- 
itatively with our data but do not show any low-spin be- 
haviour, since they do not include the region very close to 
FCC. Furthermore, those data were obtained for slightly 
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larger (experimental) volume ||6] which stabihses the high- 
moment solution. 
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Figure 5 Calculated MCA (in VCA) as a function of 
cobalt concentration x and in-plane lattice parameter a: (a) 
CVBP-geometry (Model I) and (b) EBP-geometry (Model 
II). Dark (brown and red) areas denote large MCA values 
with an easy axis along [001], light (yellow) areas denote 
small MCA. Black Unes indicate zero MCA. 

3.3 Magnetocrystalline anisotropy: chemically 
ordered films We now turn to the key quantity for mag- 
netic materials applications, the MCA. Recall, that for 
chemically ordered Fe-Co layers BNEH predicted a peak 
in Z\£'mca with a maximum height of about 800 /ieV in 
the vicinity of x = 0.6; 1.20 < c/a < 1.25. Figure 
shows our data for AEmca{x, a), evaluated with Model 
I (VCA, CVBP). They agree well with the data published 
by BNEH which were confirmed quahtatively by experi- 
ments liiiioiiiiiiiii. 

To check whether consideration of a- and x-specific re- 
laxation of the film perpendicular to the film plane might 
affect the MCA, we repeated the calculation for the EBP 
(Model II, Figure |5})). By comparison of the two panels 
of Figure |5] we find that both geometries yield qualita- 
tively the same results. A cut at a; = 0.6 (approximately 
through the peak maximum) shows, that both structure 
models yield almost the same maximum MCA energy of 
800 /ieV, see Figure |6] The main difference consists of a 
shift of the MCA-maximum from a = 2.68 A (CVBP) to 
a — 2.64 A (EBP). Further, the EBP-peak has a somewhat 
smaller width. In both cases, the c/a ratio lies between 1.22 




Figure 6 Calculated MCA (in VCA) for cobalt concentra- 
tion X = 0.6 as a function of the in-plane lattice parameter 
a: comparison between CVBP (Model I) and EBP (Model 
II) (same data as in Figure |5] but on a finer grid). 



and 1.24 (as found by BNEH). Thus, the different optimum 
a-value originates from the volume difference (experimen- 
tal vs. calculated volume). 



700 

600 

^500 
o 

^ 400 
> 

A 300 

< 
g 

Kf 200 
< 

100 




X = 0.500 

n • \ ' \ • r 



' Pd(OOl) 



X = 0.625 

-1 I r- 



X = 0.750 

I •■■•VGA 

DtiEAVG 

LI 



2.5 2.6 2.7 2.8 



2.5 2.6 2.7 2.8 2.5 2.6 2.7 2.8 

a [A] 



Figure 7 Calculated MCA for cobalt concentrations x = 
0.5, 0.625, and 0.75 as a function of the in-plane lattice pa- 
rameter a: comparison between chemically ordered struc- 
tures. Models II, III, and IV (VCA, LIq, and Ll2-derived 
structure, respectively) and ensemble-averaged supercells 
simulating disorder (EAVG, Model V). The VCA data are 
the same as in Figure|5j), but were evaluated here on a finer 
grid. The red cross denotes recent experimental data ob- 
tained forFeo.5Coo.5|Pd(001) HO). 

3.4 Magnetocrystalline anisotropy: the effect of 
chemical disorder Finally, we investigate how chemi- 
cal disorder of the films might influence their magnetic 
anisotropy. To this end, VCA results are compared with 
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results for two other structure types with perfect chemical 
order on the one hand and with an ensemble average sim- 
ulating disorder on the other hand. We recall, that VCA 
is a model imposing perfect chemical order. The calcula- 
tions were performed for the most interesting concentra- 
tion range between x = 0.5 and x = 0.75. 

Figure Q compiles calculated data for all four men- 
tioned structure models. Obviously, the specific type of 
crystallographic structure is, in the present case, of minor 
importance. What matters is chemical order This finding is 
not at all trivial, since the magnetic anisotropy is known to 
be sensitive to the detailed electronic structure and, hence, 
to the specific crystal structure. It was already mentioned 
by BNEH that LIq results for Z\_Emca closely resemble 
the VCA results for x = 0.5. We here confirm this point 
(left panel of Figure|7]i and add, that also by using an LI2- 
derived structure the VCA data {x = 0.75) are rather well 
reproduced (right panel). 

However, a striking difference is immediately visible 
between the MCA of the chemically ordered structure 
models and the MCA of Model V, the chemically disor- 
dered case. The latter yields AEmca peaks smaller than 
those of the ordered structures by factors of 1.5 to 3, 
depending on x. We suggest, that the relatively low exper- 
imental values (compared with the prediction by BNEH) 
might be caused by the lack of chemical order 

As a concluding remark, the experimental A'„ value for 
Feo.sCoo.s |Pd(001) IfTOl . red cross in the left panel of Fig- 
ure |7] seems to compare quite well with our VCA or L 1 
data. One has to bear in mind, however, that renormalisa- 
tion of the calculated curves to the experimental volume 
would shift all calculated data by 0.04 A to the right, thus 
providing coincidence of the experimental value with the 
simulated disorder. 

4 Conclusions We have shown, that neither detailed 
structural relaxation nor the specific structure type have an 
important influence on the magnetic and elastic properties 
of strained Fe-Co films. On the other hand, chemical disor- 
der can reduce the magnetic anisotropy energy by factors 
of 1.5 to 3. We suggest that this might be a reason for the 
relatively small Ku-values measured, hitherto. 

Epitaxial growth of thick, highly strained Fti^xCox 
films on substrates with lattice parameter a is predicted to 
be most stable in the parameter range 0.67 < x < 0.86; 
2.64 A < a < 2.76 A, independent of chemical order The 
maximum possible values of are expected for almost 
the same range of in-plane lattice parameters, 2.60 A < a 

< 2.75 A but somewhat lower Co concentration, 0.45 < x 

< 0.70. 

Based on these findings, we suggest two promising 
routes to fabricate strained Fe-Co films with large perpen- 
dicular anisotropy: 

(1) Preparation of relatively thin ordered Fe-Co LIq- 
films {x = 0.5) on a substrate with a « 2.64 A. The 
maximum thickness of these films is probably limited to 



about 10 ML by a relatively large strain energy. On the 
other hand, very large A'„ values up to 500 /ieV per atom 
are expected (Figure |7]i. Thus, a large volume-integrated 
anisotropy can be achieved with relatively thin films. 

(2) Preparation of relatively thick films close to the sad- 
dle point of the strain energy, e.g., x = 0.75 and 2.60 A < 
a < 2.65 A. For these parameters, Ku of chemically disor- 
dered films is limited to about 100 /ieV per atom, too small 
a value to overcome the in-plane shape anisotropy. There 
is, however, a chance to enhance the magnetic anisotropy 
of these films by careful annealing; calculations performed 
by a group including Manfred Fahnle point to the pos- 
sibility to stabilise an ordered CoaFe phase by epitaxial 
strain |231. 

We hope that our results will add another important 
route to the experimental search for Fe-Co based mag- 
netic recording media: the preparation of chemically or- 
dered strained films. 
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